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We investigate the role played by fast quenching on the decay of metastable (or false vacuum) 
states. Instead of the exponentially-slow decay rate per unit volume, Fhn ~ exp[— i?i,/fcsT] {Eb 
is the free energy of the critical bubble), predicted by Homogeneous Nucleation theory, we show 
that under fast enough quenching the decay rate is a power law Frn ~ [Et/kBT]~^ , where B is 
weakly sensitive to the temperature. For a range of parameters, large-amplitude oscillations about 
the metastable state trigger the resonant emergence of coherent subcritical configurations. Decay 
mechanisms for different Eb are proposed and illustrated in a (2-|-l)-dimensional scalar field model. 

PACS numbers: ll.lO.Lm, 64.60.Qb, 05.45.Xt, 98.80.Cq 



Introduction. Few topics in physics have the range 
of applicability of first order (or discontinuous) phase 
transitions. From materials science to particle physics 
and cosmology, the fact that a large number of physi- 
cal systems can be described as having two phases sep- 
arated by an energy barrier has been an active topic 
of research for decades Q |^ [H. Much of the theo- 
retical work in this field is derived in one way or an- 
other from the theory of Homogeneous Nucleation (HN) 
Q 01 . HN assumes that the system is initially local- 
ized in a spatially-homogeneous metastable state, that is, 
that only small fiuctuations about the local equilibrium 
state, 0m, exist: one computes the partition function 
Z = —TlnJ^ summing only over quadratic fluctuations 
about (f>m- is the free energy, given by the path integral 
/ X'0exp[— i?[(/)]/r], where E[(j)] is the free energy func- 
tional of configuration (f). [ks = c = h = 1 throughout.] 
In relativistic field theory, false vacuum decay has been 
examined both at zero [3and finite temperature 6] by a 
large number of authors [7l . In general, the HN approxi- 
mation is adopted from the start. At finite temperatures, 
one uses the well-known exponential decay rate per unit 
volume, Fhn ^ T^'^+i) exp[-i;f,/r], where Et is the free 
energy barrier for the decay, or the energy of the critical 
bubble or bounce, 0b (r), the solution to the equation 
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with appropriate boundary conditions. A prime denotes 
derivative with respect to the d-dimensional radial co- 
ordinate (sphericity is energetically favored) and V[(j)] is 
the effective potential that sums over thermal and quan- 
tum contributions when applicable. At T = 0, one has a 
purely quantum vacuum decay, and the pre-factor T^'^^^'> 
is roughly approximated by M^'^+^\ the relevant mass 
scale, while Eb/T is substituted by SE[(f>b], the {d + 1)- 
dimensional Euclidean action of the bounce configura- 
tion. 



In the present work we examine what happens if one 
relaxes the HN approximation that the initial state is 
well-localized about equilibrium. We subject the sys- 
tem to an instantaneous quench, equivalent to a sudden 
change of potential from a single to an asymmetric double 
well. This should be contrasted with the work of ref. Q 
which studies quenches in models without a barrier sepa- 
rating symmetric and broken-symmetric states, and thus 
with spinodal decomposition dynamics. Although in this 
first study we will only consider instantaneous quenches, 
we expect our results to carry on at least partially to 
slower quenches, so long as the quenching rate Tquonch is 
faster than the relaxation rate of the field's zero mode, 
tq. Why the field's zero mode? With the longest wave- 
length it is the slowest to equilibrate: as Tquench tq, 
the system will remain in equilibrium. For appropriate 
choices of parameters, the rapid quench will induce large- 
amplitude oscillations of the field's zero mode |^. Due to 
the nonlinear potential, energy will be transferred from 
the zero mode to higher fc-modes. As observed in refer- 
ence , this transfer of energy results in the synchronous 
emergence of oscillon-like configurations 1^. [We urge 
the reader to consult reference 9] for details.] For small 
enough double-well asymmetry, these localized field con- 
figurations act as precursors for the nucleation of a criti- 
cal bubble, greatly reducing the decay time-scale. In the 
simulations we examined, the critical bubble emerges as 
two or more subcritical oscillons coalesce, or, for larger 
asymmetries, as a single oscillon becomes critically un- 
stable to growth. If the asymmetry is too large, the field 
crosses directly to the global minimum. 
The Model. Consider a (2-|-l)-dimensional real scalar 
field (or scalar order parameter) i)f)(x, t) evolving under 
the infiuence of a potential V ((/>). The continuum Hamil- 
tonian is conserved and the total energy of a given field 
configuration 0(x, t) is. 
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where V{(j)) = ^4>^ - ^4>^ + jcj)^ is the potential energy 
density. The parameters m, a, and A are positive defi- 
nite and temperature independent. It is helpful to intro- 
duce the dimensionless variables 4>' = (j) VA/to, x' = xm, 
t' = tm, and a' = a/ (m^/X) (We will henceforth drop the 
primes). Prior to the quech, a = and the potential is 
an anharmonic single well symmetric about (f> — 0. The 
field is in thermal equilibrium with a temperature T. At 
the temperatures considered, the fluctuations of the field 
are well approximated by a gaussian distribution, with 
((jp ) = aT (a = 0.51 and can be computed numerically 
[llj). As such, within the context of the Hartree approx- 
imation |l2l |. the momentum and field modes in /c-space 
can be obtained from a harmonic effective potential, and 
satisfy (|7f(fc)H - T and (|0(k)|2) = ^5^, respec- 
tively. The Hartree mass rnj^ = 1 + 1(0^) depends on 
the magnitude of the fluctuations (and thus T). Here- 
after we will refer to a particular system by its initial 
temperature. All results are ensemble averages over 100 
simulations. 

If a 7^ 0, the Z2 symmetry is explicitly broken. When 
a ~ 1.5 = ac, the potential is a symmetric double- well 
(SDW), with two degenerate minima. The quench is im- 
plemented at t = by setting a > ac, whereby the poten- 
tial is asymmetric (ADW) about the barrier separating 
the two minima. The Hartree approximation gives an ac- 
curate description of the evolution of the area-averaged 
field (/)avo(i) and its fluctuations for early times after the 
quench, in which the distribution of fluctuations remains 
gaussian and the dynamics are governed by an effective 
potential, 

- f Cc + ^'/>L- (3) 

The quench shifts the local minimum from 0avc = to a 
positive value, but also introduces a new global minumum 
to the system. [See inset in Fig. |21] 
Nucleation of Oscillons. The quench induces oscilla- 
tions in 0avo about the new local minimum, which even- 
tually dampen due to nonlinear scattering with higher k- 
modes. At early times small fluctuations satisfy a Math- 
ieu equation in /c-space 

'^>=-[fc'+Kff['/'avc(t)]]'50, (4) 

and, depending on the wave number and parametric os- 
cillations of 4>a,Me(t), can undergo exponential amplifica- 
tion e^*). In Fig. we show the lines of constant 
amplification rate for different wave numbers and tem- 
peratures when a = Qfc and Voft is defined by the initial 
thermal distribution. At low temperatures T < 0.13, 
no modes are ever amplified. As the temperature is in- 
creased, so is the amplitude and period of oscillation in 
0ave, gradually causing the band < fc < 0.48 to resonate 
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FIG. 1: Lines of constant amplification rate 77 for small- 
amplitude modes at various temperatures, beginning with 
T?min = 2.8 X 10~^ for the bottom-most contour and increasing 
in increments of At; = 1.3 x 10"'^. 

and grow. A full description of the coupled dynamics of 
4'a.vc(t) and d<p{:x.,t) when a = ac is given in ref. Ill|. 
Furthermore, for large enough temperatures (T > 0.13) 
large-amplitude fiuctuations about the zero mode probe 
into unstable regions where V^g < 0, which also pro- 
mote their growth. Note that this is very distinct from 
spinodal decomposition, where competing domains of the 
metastable and stable phases coarsen Ij. Instead, for the 
values of T and a considered, (/)avc continues to oscillate 
about the metastable minimum until a critical bubble of 
the stable phase grows to complete the transition. 

These two processes result in the synchronous emer- 
gence of oscillon-like configurations (Fig. 3 in ref. Q), 
long-lived time-dependent localized field configurations 
which are well-described by gaussian profiles, (/)osc(^, f) ~ 
(j>a{t) exp[—r'^ / R"^] To strengthen our argument, 

note that within this gaussian ansatz, an oscillon is com- 
prised by modes within the band < k < 2/ R. One 
of us has recently shown that, in d dimensions and for 
a potential V, the radius of an oscillon satisfies R^ > 
d/[i(23/V3)'^(y"')VV^^^ - V"] 13]. For the potential 
of eq. 121 and d — 2, we obtain that the related band of 
wave numbers is, < /c < 0.66. Referring back to Fig. ^ 
the reader can verify that these are also approximately 
the modes excited by parametric resonance. 
Resonant Nucleation. Having established that oscil- 
lons emerge after the quench, we can examine their role 
as precursors of metastable decay. Unless the potential 
has a large asymmetry, oscillons are typically sub-critical 
fiuctuations; as will be discussed below, a critical nucleus 
may appear only due to the coalescence of two or more 
oscillons. It should be clear, however, that their appear- 
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FIG. 2: The evolution of the order parameter (/f>avc(i) at 
T = 0.22 for several values of the asymmetry. From left to 
right, a = 1.746, 1.56, 1.542, 1.53, 1.524, 1.521, 1.518. The in- 
set shows Veff for the same values. 
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FIG. 3: Decay time-scale trn as a function of critical nucle- 
ation effective free-energy barrier Ei,/T at T = 0.18, 0.20, and 
T = 0.22. The best fits (dashed lines) are power-laws with 
exponents B ~ 3.762, 3.074, and 2.637, respectively. 



ance renders the homogeneity assumption of HN theory 
inapplicable: the metastable background is far from ho- 
mogeneous and the critical energy barrier must be rcnor- 
malized [l5| . In other words, a rapid quench or cooling 
leads to departures from the usual HN assumptions. As 
we describe next, the decay rate of the quenched system 
may be much faster than what is predicted by HN theory. 

In Fig. 121 we show the evolution of the order pa- 
rameter 0avo(^) a-s a function of time for several val- 
ues of asymmetry, 1.518 < a < 1.746, for T = 0.22. 
Not surprisingly, as a — > ac = 1.5, the field remains 
longer in the metastable state, since the nucleation en- 
ergy barrier i?f, —> oo at etc. However, a quick glance at 
the time axis shows the fast decay time-scale, of order 
10"'^"^. For comparison, for 1.518 < a < 1.56, homo- 
geneous nucleation would predict nucleation time-scales 
of order ~ lO^^ > trn ~ exp[£;f,/r] > lO^^ (in di- 
mensionless units). [The related nucleation barriers with 
the effective potential are E}j{a — 1.518) — 14.10 and 
£"{,(« = 1.56) = 5.74.] While for smaller asymmetries 
4>ave{t) displays similar oscillatory behavior to the SDW 
case before transitioning to the global minimum, as a is 
increased the number of oscillations decreases. For large 
asymmetries, a > 1.746, the entire field crosses over to 
the global minimum without any nucleation event, re- 
sulting in oscillations about the global minimum. The 
inset in Fig. |21 shows that the barrier is just low enough 
for this to occur. 

In Fig. 13 we show the ensemble-averaged nucleation 
time-scales for resonant nucleation, trn, as a function of 
the nucleation barrier (computed with eq. EJ, Eb/T, for 
the temperatures T = 0.18, 0.20, and 0.22. [For tem- 
peratures above T = 0.26 one is in the vicinity of the 



critical point in which no barrier exists.] The nucleation 
time was measured when (^ave crosses the maximum of 
Vctf. The best fit is a power law, trn oc (Eh/T)^ , with 
B = 3.762 ± 0.016 for T = 0.18, B = 3.074 ± 0.015 for 
T = 0.20, and B = 2.637 ±0.018 for T = 0.22. This sim- 
ple power law holds for the same range of temperatures 
where we have observed the synchronous emergence of os- 
cillons. It is not surprising that the exponent B increases 
with decreasing T, since the synchronous emergence of 
oscillons becomes less pronounced and eventually van- 
ishes. In these cases we should expect a smooth transi- 
tion into the exponential time-scales of HN. We present 
below what we believe is the mechanism by which the 
transition completes for different nucleation barriers. 

First, for a — > ac, the radius of the nucleation bub- 
ble diverges, ^ oo. When fast quenching induces 
large-amplitude fluctuations of the field's zero mode, the 
system doesn't approach the global minimum through a 
random search in configuration space as is the case in 
HN. Instead, we argue that oscillons will induce the nu- 
cleation of a critical fluctuation. The way in which this 
happens depends on the magnitude of the nucleation bar- 
rier: for nearly degenerate potentials, Uc < a < ai, the 
critical nucleus has a much larger radius than a typical 
oscillon; it will appear as two or more oscillons coalesce. 
We call this Region I, defined for Ri, > 2i?osc, where Rns r. 
is the minimum oscillon radius computed from ref. [iSj. 
Fig. 0] illustrates this mechanism. Two oscillons, labeled 
A and B, join to become a critical nucleus. They dif- 
fuse through the lattice and form bound states, some- 
what as in kink-antikink breathers in Id field theory 
[l^ . [The interested reader can see simulation movies 
at Mtp://www.dartmouth.edu/^cosmos/oscillons ] We 





FIG. 4: Two oscillons coalesce to form a critical bubble. First 
two frames from top show oscillons A and B. Third and fourth 
frames shows A and B coalescing into a critical bubble. Fi- 
nal frame shows growth of bubble expanding into metastable 
state. 

are currently attempting to estimate the diffusion and 
coalescence rate of oscillons on the lattice so that we can 
compute the power law decay rate analytically. 

As a is increased further, the radius of the critical nu- 
cleus decreases, approaching that of an oscillon. In this 
case, a single oscillon grows unstable to become the criti- 
cal nucleus and promote the fast decay of the metastable 
state: there is no coalescence. We call this Region II, 
< Q^iIj ^6 < 2i?osc- This explains the small num- 
ber of oscillations on 0ave(O is increased [cf. Fig. 
12]. To corroborate our argument, in Fig. |Slwe contrast 
the critical nucleation radius with that of oscillons as ob- 
tained in ref. different values of effective energy 
barrier and related values of a at T = 0.22. The crit- 



FIG. 5: Radius of critical bubble (Rt) and twice the minimum 
oscillon radius (2i?osc) as a function of its energy barrier and 
related values of a at T = 0.22. For a > 1.547 one cannot 
easily distinguish between an oscillon and a critical bubble. 



ical nucleus radius Rh is equal to 2i?osc for a — 1.547. 
This defines the boundary between Regions I and II: for 

^ ai SL single oscillon may grow into a critical bubble. 
Finally, for a > au = 1.746 the field crosses over to the 
global minimum without any nucleation event. 

How will the efficiency of the mechanism decrease as 
Tquonch —^tq? What happens when the quench is induced 
by cooling as, for example, in the early universe? Does 
the power law behavior obtained here still hold for d = 31 
We intend to address these and related questions in the 
near future. 
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